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' The limitations of the recently proposed new method of numerical modelling of 

^ l' Bose-Einstein correlations (BEC) are explicitly demonstrated. It is then argued 

that BEC should still be considered as emerging from the correlations of fluctua- 
tions, however they have to be modelled first in any Monte Carlo event generator 
(MCEG) and not added a posteriori to the existing output of some MCEG. 



Oh! 
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. Recently we have proposed new method of numerical modelling of BEC 

$—( ' which makes use of the fact that BEC can be regarded as correlations of fluc- 

^ tuations Our aim was to provide a fast algorithm, which could be used 

on the even-by-event basis together with the output of other MCEG's and 
which would introduce a characteristic structure of the correlation function 
C2{Q) between them ^. The main idea was ^'^ to resign from the initial 
allocation of charges to the particles provided by a given MCEG and to 
perform then the new allocation in a way which would result in a desired 
bunching of like particles in the phase space**. As a result one gets a number 
of (what we have called ^) elementary emitting cells (EEC), Nceii, each of 
them containing particles, Upart, of a given charge distributed according to 
geometrical distribution (in order to model the bosonic character of pro- 
duced secondaries '^''^). Those cells (and parameters connected with their 
formation) represent therefore the main object of interest in our approach. 



^At the same time all positions of the initial particles in the phase-space, as well as the 
number of particles of given charge in an event, remain intact l'^'^. 
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It means then that their number and the (mean) multipUcity of secondaries 
they contain are the main factors dictated the shape of C^iQ) (and there- 
fore also the 'size' parameter R and 'chaoticity' parameter A usually used 
to fit C2{Q) ^'^)- In Fig. 1 we show that it is, in principle, possible to fit 
experimental data. Our reservation comes from the fact that so far only 
one sample of e+e" data were fitted by only two types of MCEG (cf. ^ for 
details) and using only direct pions. But it shows also explicitly the price 
to be paid for this, namely the anticorrelations between the unlike-particles 
showing up^. 
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Figure 1. Comparison with DELPHI 
data presented in shown together with 
the results for the unlilce-pairs for CAS 
model ^. 



Figure 2. Comparison with DELPHI 
data presented in ^ are shown together 
with the results for the unlike-pairs for 
JETSET model 3. 



Such anticorrclation feature is, however, after a bit of thinking, to be 
expected. The reason is that we start, in a given event, with a priori uncor- 
related particles provided by a MCEG. Our algorithm shifts now charges in 
a fixed pool of particles in this event. It is constructed in such a way that the 
EEC of different charges essentially do not overlap in the phase space (al- 
though, in principle there is nothing preventing them to do so; but we have 
checked that allowing for substantial overlap in order to diminish the anti- 
correlation effect weakness dramatically the obtained BEC signal as well). 
The possible way out is the reparametrization of the original MCEC in such 
way as to have positive charge-particles correlations from the very begin- 
ning which are then washed out by the action of our algorithm. Another, 



''We are deeply indebted to Prof. W.Kittel for suggesting to us this simple check before 
proceeding with any further, more involved, development of our algorithm. 
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Figure 3. (a, 6) Application of our toy model with two types of weights P (defining 
the formation of the EEC's) to fit DELPHI data on e+e^ annihilation. Panels (c) 
and (d) show the corresponding distributions of EEC, PiN^g^u), and distribution of the 
like-charge particles allocated to such cell, P{npart), with mean values and dispersions 
as indicated (corresponding also to parameters of the best fits, respectively, poissonian 
and geometrical, to P{Nceii) and P{npart), which are not shown here). Notice the small 
values of multiplicities in EEC's, (n) ~ 1.1 only. 

more physical (in our opinion) possibility, which we shall demonstrate here, 
is to construct MCEG satisfying the requirements of BEG from the very 
beginning (i.e., essentially to continue the line of reasoning proposed some 
time ago in ^). It means that, contrary to the situation in Figs. 1 and 2, 
where the original outputs of MCEG were not showing any BEG and our 
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algorithm was not changing the main characteristics of theses outputs (hke, 
for example, multiplicity distributions), one has to change from the very be- 
ginning the inclusive distributions to accommodate the characteristic BEC 
effects in them, see also discussion in Fig. 3 shows results of our first 
attempts in this direction ^ where preliminary comparison of a kind of toy 
model (presented in ^) with data at = 91.2 GeV has been shown for two 
types of the specific weight variables, P ^. Particles were selected from the 
energy W with given proportion of (+/ — /O) and allocated to EEC's in such 
a way that, as evident in Fig. 3, distribution of cells, P{Nceii), is poissonian 
whereas distribution of particles (of the same sign) in a cell, P{npart), is 
geometrical. It means then that the multiplicity distribution of all particles 
follow the Negative Binomial (NB) form with EEC playing a role of clans 
^. It is very important to realize that in order to get the observed structure 
of C2{Q) one has to allow for some uncertainties in energies of particles 
allocated to a given EEC (here given by a gaussian form with width a as 
indicated in Fig. 3, they correspond to the sizes of EEC, 6y, introduced in 
although they are not identical with them). As one can see the result of 
our "toy model" shown in Fig. 3 is quite promising and can form a basis 
of more detailed investigations (which we plan to pursue*^). 
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°In particular one should improve the approximate for a moment (although kept on the 
same numerical level as in conservation of charges and energy- momenta in a given 
event. 



